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Abstract 

A twistor construction of the hierarchy associated with the hyper-Kahler equations on a metric (the 
anti-self-dual Einstein vacuum equations, ASDVE, in four dimensions) is given. The recursion operator 
R is constructed and used to build an infinite-dimensional symmetry algebra and in particular higher 
flows for the hyper-Kahler equations. It is shown that R acts on the twistor data by multiplication with 
a rational function. The structures are illustrated by the example of the Sparling- Tod (Eguchi-Hansen) 
solution. 

An extended space-time Af is constructed whose extra dimensions correspond to higher flows of the 
hierarchy. It is shown that A/" is a moduli space of rational curves with normal bundle 0{n) © 0{n) in 
twistor space and is canonically equipped with a Lax distribution for ASDVE hierarchies. The space Af 
is shown to be foliated by four dimensional hyper-Kahler slices. 

The Lagrangian, Hamiltonian and bi-Hamiltonian formulations of the ASDVE in the form of the 
heavenly equations are given. The symplectic form on the moduli space of solutions to heavenly equations 
is derived, and is shown to be compatible with the recursion operator. 

1 Introduction 

Roger Penrose's twistor theory gives rise to correspondences between solutions to differential equations on 
the one hand and unconstrained holomorphic geometry on the other. The two most prominent systems of 
nonlinear equations which admit such correspondences are the anti-self-dual vacuum Einstein equations 
(ASDVE) which in Euclidean signature determine hyper-Kahler metrics, and the anti-self-dual Yang- 
Mills equations (ASDYM) Richard Ward observed that many lower-dimensional integrable 
systems are symmetry reductions of ASDYM. This has led to an overview of the theory of integrable 
systems [Q, which provides a classification of those lower-dimensional integrable systems that arise as 
reductions of the ASDYM equations and a unification of the theory of such integrable equations as 
symmetry reduced versions of the corresponding theory of the ASDYM equations. In [Q, Lagrangian 
and Hamiltonian frameworks for ASDYM were described together with a recursion operator. This leads 
to the corresponding structures for symmetry reductions of the ASDYM equations. 

In this paper we investigate these structures for the second important system of equations — the ASDVE 
or hyper-Kahler equation (this system also admits known integrable systems as symmetry reductions ) . 
We shall give a twistor-geometric construction of the hierarchies associated to the ASDVE in the 'heavenly' 
forms due to Plebahski p5| . In this context it is more natural to work with complex (holomorphic) metrics 
on complexified space-times and so we use the term ASDVE equations rather than hyper-Kahler equations. 
Our considerations will generally be local in space-time which will be understood to be a region in C^. 

In Section || we summarise the twistor correspondences for flat and curved spaces. We establish a 
spinor notation (which will not be essential for the subsequent sections) and recall basic facts about the 
ASD conformal condition and the geometry of the spin bundle. In Section || the recursion operator R 
for the ASDVE is constructed as an integro- differential operator mapping solutions to the linearised 
heavenly equations to other solutions. We then use this to give an alternate development of the twistor 
correspondence by using R to build a family of foliations by twistor surfaces. We show that R corresponds 
to multiplication of the twistor data by a given twistor function. We then analyse the hidden symmetry 
algebra of the ASDVE, and use the recursion operator to construct Killing spinors. We illustrate the 
ideas using the example of the Sparling-Tod solution and show how R can be used to construct rational 
curves with normal bundle 0(1) © 0{1) in the associated twistor space. 

In Section ^ we give the twistor construction for the ASDVE hierarchies. The higher commuting flows 
can be thought of as coordinates on an extended space-time. This extended space-time has a twistor 
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correspondence: it is the moduli space of rational curves with normal bundle 0{n) 0{n) in a twistor 
space. This moduli space is canonically equipped with the Lax distribution for ASDVE hierarchies, and 
conversely that truncated hierarchies admit a Lax distribution that gives rise to such a twistor space. The 
Lax distribution can be interpreted as a connecting map in a long exact sequence of sheaves. In Section ^ 
we investigate the Lagrangian and Hamiltonian formulations of heavenly equations. The symplectic form 
on the moduli space of solutions to heavenly equations will be derived, and is shown to be compatible 
with the recursion operator. 

We end this introduction with some bibliographical remarks. Significant progress towards understand- 
ing the symmetry structure of the heavenly equations was achieved by Boyer and Plebahski Q who 
obtained an infinite number of conservation laws for the ASDVE equations and established some connec- 
tions with the nonlinear graviton construction. Their results were later extended in papers of Strachan 
and Takasaki ^ . The present work is an extended version of P, |9| p^. 

2 Preliminaries 
2.1 Spinor notation 

We work in the holomorphic category with complexified space-times: thus space-time is a complex 
four-manifold equipped with a holomorphic metric g and compatible volume form j/. 

In four complex dimensions orthogonal transformations decompose into products of ASD and SD 
rotations ^ 

50(4, C) = iSL{2, C) X SL{2, C))/Z2. (2.1) 

The spinor calculus in four dimensions is based on this isomorphism. We use the conventions of Penrose 
and Rindler [Q. Indices will generally be assumed to be concrete unless stated otherwise: a,b,..., 
a — 0, 1 ... 3 are four-dimensional space-time indices and A, B, . . . , A' , B' , . . ., A — 0, 1 etc. are two- 
dimensional spinor indices. The tangent space at each point of A4 is isomorphic to a tensor product of 
the two spin spaces 

T''M^S'^(g>S'^'. (2.2) 

The complex Lorentz transformation — > A"hV^'', A'^^A'^dgac — gbd, is equivalent to the composition 
of the SD and the ASD rotation 

where X^b and A'* b' are elements of iS'L(2,C) and iS'L(2,C). 

Spin dyads (o"^,i"*) and {o^ ) span and respectively. The spin spaces S"^ and are 
equipped with symplectic forms eab and ea'B' such that Eqi — £o'1' = 1- These anti-symmetric objects 
are used to raise and lower the spinor indices. We shall use normalised spin frames so that 

O^^C- _ ^B^C ^ ^BC^ ^B'^C „ ^B'^C ^ ^B'C'^ 

Let e^'^ be a null tetrad of 1-forms on M and let V aa' be the frame of dual vector fields. The orientation 
is given by fixing the volume form 

Apart from orientability, M must satisfy some other topological restrictions for the global spinor fields 
to exist. We shall not take them into account as we work locally in M. 

The local basis and E^'-^' of spaces of ASD and SD two-forms are defined by 

^AA' ^ ^BB' ^ ^AB^A'B' ^ ^A' B' ^AB _ (2.3) 

The first Cartan structure equations are 

de^-^' =e^^' AL^B + e^^' AT-^'s,, 

where T ab and T a'B' are the SL(2, C) and SL(2, C) spin connection one-forms. They are symmetric in 
their indices, and 

= ^CCABS^^ , ^A'B' — ^CC'A'B'S^^ , ^CC'A'B' — OA'^CC-B' — iA'^ CC'OB' ■ 



The curvature of the spin connection 



R^B = dV^B + r^c A T^B 



decomposes as 



and similarly for b' ■ Here R is the Ricci scalar, '^abA'B' is the trace-free part of the Ricci tensor Rab, 
and Cabcd is the ASD part of the Weyl tensor 

Cabcd — £A'B'£C'D'CaBCD + £ AB^CdC A' B' C D' ■ 



2.2 The flat twistor correspondence 

The flat twistor correspondence is a correspondence between points in complexified Minkowski space, 
(or its conformal compactification) and holomorphic lines in CP"^. 

The flat twistor correspondence has an invariant formulation in terms of spinors. A point in has 



position vector with coordinates {w,z,x,y). The isomorphism (2.2) is realised by 



x^^' := y ^zj' so that g = eabsa' B'dx'^'^' dx^^' . 

A two-plane in is null if g{X, Y) —Q for every pair (X, Y) of vectors tangent to it. The null planes can 
be self-dual (SD) or anti self-dual (ASD), depending on whether the tangent bi- vector X f\Y is SD or ASD. 
The SD null planes are called a-planes. The a-planes passing through a point in are parametrised by 
A = tto'/tti' S CP^. Tangents to a-planes are spanned by two vectors 

^.-"^ (2.4) 

which form the kernel of tta'T^B''^^ ^ The set of all a-planes is called a projective twistor space and 
denoted VT. For it is a three-dimensional complex manifold bi- holomorphic to CP'^ — CP^ . 

The five complex dimensional correspondence space := C* x CP^ fibres over C''^ by (a;^'^ , A) x^'^ 
and over VT with fibres spanned by La- Twistor functions (functions on PT) pull back to functions on 
T which are constant on a-planes, or equivalently satisfy La] — 0. 

Twistor space can be covered by two coordinate patches U and ?7, where J7 is a complement oi\ — co 
and [/ is a compliment of A = 0. If (/i°, A) are coordinates on U and {Jp , P'^jX) are coordinates on U 
then on the overlap 

^" = //A, /i^ = ^i^/X, A = 1/A. 
The local coordinates {ijP, fi^, A) on VT pulled back to T are 

fi^ — w + Xy, /i^ = z — Xx, X. (2-5) 

We can introduce homogeneous coordinates on the twistor space 

{uj^,TTA') = (w",a;\7ro',7ri') (^°7ri/, ^^tti', Atti/, tti/). 

The point x"^^ G C lies on the a-plane corresponding to the twistor (lu^, tta') G VT iff 

Lo^ ^x'^^'tta'. (2.6) 

For TTA' ^ and {uj"^, tt^') fixed. The solution to ( |2.6| ) is a com plex two plane with tangent vectors of the 
form TT^ for all . Alternatively, if we fix x^^' , then ( ^.6| ) defines a rational curve, CP"'^, in VT with 
normal bundle C'(1)®0(I)|^ Kodaira theory guarantees that the family of such rational curves in VT is 
four complex dimensional. There is a canonical (quadratic) conformal structure ds^ on C*: the points p 
and q are null separated with respect to ds^ in iff the corresponding rational curves Ip and Iq intersect 
in VT at one point. 



2.3 Curved twistor spaces and the geometry of the primed spin bundle. 

Given a complex four-dimensional manifold AA with curved metric g, a twistor in M. is an a-surface, 
i.e. a null two-dimensional surface whose tangent space at each point is an a plane. There are Frobenius 
integrability conditions for the existence of such a-surfaces through each a-plane element at each point and 
these are equivalent, after some calculation, to the vanishing of the self-dual part of the Weyl curvature, 
Ca'B'C'D'- Thus, given Ca'B'C'D' — 0, we can define a twistor space VT to be the three complex 
dimensional manifold of a-surfaces in A1 . If <? is also Ricci flat then VT has further structures which are 
listed in the Nonlinear Graviton Theorem: 

^ Here 0{;n) denotes the line bundle over CP^ with transition functions A~" from the set A 7^ oo to A 7^ (i.e. Chern class 

n). 



Theorem 2.1 (Penrose |23|) There is a 1-1 correspondence between complex ASD vacuum metrics on 
complex four-manifolds and three dimensional complex manifolds VT such that 

• There exists a holomorphic projection ^ : VT — > CP^ 

• VT is equipped with a four complex parameter family of sections of /i each with a normal bundle 
0{\) © 0{\), (this will follow from the existence of one such curve by Kodaira theory), 

• Each fibre of ^ has a symplectic structure T,x G T{A^{ii~^{\)) (E) 0(2)), where A € CP"'^. 

To obtain real metrics on a real 4-manifold, we can require further that the twistor space admit an 
anti-holomorphic involution. 

The correspondence space T — M. X CP^ is coordinatized by (x, A), where x denotes the coordinates 
on M and A is the coordinate on CP^ that parametrises the a-surfaces through x \n M. We represent 
T as the quotient of the primed-spin bundle with fibre coordinates tt^' by the Euler vector field 
T ~ TT^ /diT^ . We relate the fibre coordinates to A by A = -kq' /ttv . A form with values in the line bundle 
0{n) on T can be represented by a homogeneous form a on the non-pro jective spin bundle satisfying 

T J a = , Cya = na. 

The space T possesses a natural two dimensional distribution called the twistor distribution, or Lax 
pair, to emphasise the analogy with integrable systems. The Lax pair on arises as the image under the 
projection TS^ — s- TT of the distribution spanned by 



on TS^ where the Oaa' are a null tetrad for the metric on Al, and Taa'B'C are the components of the 
spin connection in the associated spin frame {Daa' + ^ aa' b'C'T^^ aFT horizontal distribution on 

Sa')- We can also represent the Lax pair on the projective spin bundle by ^ 

LA^{TTy^){7:^'dAA'+ fAdx), where /a = (^i^')rAA'S'C'^'^V^V'^'. (2.8) 

The integrability of the twistor distribution is equivalent to Ca'B'C'D' = 0, the vanishing of the self-dual 
Weyl spinor. When the Ricci tensor vanishes also, a covariant constant primed spin frame can be found 
so that Taa'B'C = 0. We assume this from now on. 

The projective twistor space VT arises as a quotient of J- by the twistor distribution. With the Ricci 
flat condition, the coordinate A descends to twistor space and tta' descends to the non-projective twistor 
space. It can be covered by two sets, U = {|A| < I -I- e} and U = {|A| > 1 — e}. On the non-projective 
space we can introduce extra coordinates lo"^ of homogeneity degree one so that (w'*, tt^'), tta' are 
homogeneous coordinates on U and similarly (a)"*, tt^i'), tt^/ ^ oa') on JJ . The twistor space VT is then 
determined by the transition function 

The correspondence space has the alternate definition 

T = VTy. M\zei. = M x CF^ 

where Ix is the line in VT that corresponds to x E M and Z E VT lies on Ix- This leads to a double 
fibration 

T -^VT. (2.9) 

The existence of La can also be deduced directly from the correspondence. From |2^, points in M 
correspond to rational curves in VT with normal bundle 0^(I):=0(1)®0(1). The normal bundle to l^ 
consists of vectors tangent to x (horizontally lifted to T^x,\)^) modulo the twistor distribution. Therefore 
we have a sequence of sheaves over CP^ 

— > D — > — > 0-^(1) — > 0. 

The map C"* — > 0^(1) is given by — > y-^^ tt^, . Its kernel consists of vectors of the form tt"* A"* 
with A"* varying. The twistor distribution is therefore D = 0(— 1) ® S"^ and so there is a canonical 
La e T{D (g) 0(1) (g) Sa), as given in ^ 



on the projective spin bundle. In a projection 5*^ — > T we shall use the replacement formula 



^ Various powers of tti/ in formulae like (2.5) guarantee the correct homogeneity. We usually shall omit them when working 

- we shall us 

— — (2.7) 



This is because (on functions of A) 



a' a' a' 

TTq' \ TVifO — TTQ/t TV 



2.4 Some formulations of the ASD vacuum condition 



The ASD vacuum conditions Ca'B'CD' — 0, ^aba'B' = = R imply the existence of a normahsed, 
covariantly constant frame {o^ , ) of , so that Taa'B'c — 0. One can further choose an unprimed 
spin frame so that the Lax pair (2.8) consists of volume-preserving vector fields on Ai: 



Proposition 2.2 (Mason &: Newman [|lq| .) Let Vaa' = (Vqo', Vqi'Vio'Vh') be four independent 
holomorphic vector fields on a four- dimensional complex manifold M and let v he a nonzero holomorphic 
four-form. Put 



Lo 



'00' 



AVoi 



Li 



' 10' 



AVii 



Suppose that for every X G 



0, 



Here Cy denotes the Lie derivative. Then 

Oaa' = f^^'^ AA' , where 



■■= i^(Voo',Voi'Vio'ViiO, 



(2.10) 
(2.11) 



is a null-tetrad for an ASD vacuum metric. Every such metric locally arises in this way. 

In the last proposition is generalised to the hyper-Hermitian case. A choice of unprimed spin frame 
with = 1 is always possible and we shall assume this here-on so that Waa' — Vaa'- For easy reference 
we rewrite the field equations (2.11) in full 

[Vao',Vbo'] = 0, 

[VacVbi-] + [Vai',Vbo' 
[VAi',Vijr] =0. 

Let E"^ ^ be the usual basis of SD two-forms. On the correspondence space, define 



0, 



S(A) := S^'^'^A'^B' 



(2.12) 

(2.13) 
(2.14) 

(2.15) 



The formulation of the ASDVE condition dual to (2.11) is: 



Proposition 2.3 (Plebanski [25], Gindikin |12|) If a two-form of the form 

E(A) S^'^'^A'^B' 



on the correspondence space satisfies 

dhS(A) =. 0, 



I](A) A E(A) = 



(2.16) 



'-elated to E^'-^' 



where d^ is the exterior derivative holding tta' constant, then there exist one-forms 
by equation ( [g.^ j which give an ASD vacuum tetrad. 

Note that the simplicity condition in (2.16) arises from the condition that T,^ ^ comes from a tetrad. 

To construct Gindikin's two-form starting from the twistor space, one can pull back the fibrewise 
complex symplectic structure on VT — > CP^ to the projective spin bundle and fix the ambiguity by 
requiring that it annihilates vectors tangent to the fibres. The resulting two-form is 0(2) valued. (To 
obtain Gindikin's two-form one should divide it by a constant section of 0{2).) 
Put SO'o' = -5, SO'i' 



Ll> A LU 



= 2a A (5 



The second equation in ( 2.16| ) becomes 
~2z^, aAuj = aAu! = aAa = aAa — 0. 



Equations ( |2.16| ) can be seen to arise from ( p.ll| ) by observing that S(A) can be defined by 

EAB^W = v{La, Lb, ...). 



Note also that La spans a two-dimensional distribution annihilating S(A). 

by definition annihilate the twistor distribution. Define (1,1) tensors 



aB' 
^A' 



The two one- forms e^ := tta'C^^ 



^AB' 



\7aa' so that 



La 



do + X{d ~d)~ \^d2 



where (9q/ , 9° , 9q, , 9^/ ) — {8,80,82,8). If the field equations are satisfied then the Euclidean slice of 
M is equipped with three integrable complex structures given hy Ji := {1(82 — do) , {8 — 8), {82 80)} 
and three symplectic structures Ui = {{i{a — a) , ilu , {a + a)} compatible with the Ji. It is therefore a 
hyp er-K abler manifold. 



2.5 The ASD condition and heavenly equations 

Part of the residual gauge freedom in ( ^.ll| ) is fixed by selecting one of Plebahski's null coordinate 



systems. 



1. Equations (|2.13| ) and ( |2.14| ) imply the existence of a coordinate system 

{w, z, w, z) =: (w^, w^) 
and a complex- valued function Q such that 

d^n d 



dAA' = 

Equation ( |2.12| ) yields the first heavenly equation 



= (- 



d 



Kdw^dw^ dwB dw^ 



^wz^zw ^ww^zz — 1 or 



1 d'^Q 



The dual tetrad is 



-,A1' 



2 dwAdwB dw^dw 



1. 



(2.17) 



(2.18) 



— dwB 



dwAdwB 

with the flat solution VL = w^wa- The only nontrivial part of S'^ 
Kahler scalar. The Lax pair for the first heavenly equation is 



Li 



(2.19) 

is E*^ ^ = ddfl so that is a 
(2.20) 



Equations = -^I'l' = have solutions provided that satisfies the first heavenly equation 

( 2.fg ). Here ^ is a function on !F. 

2. Alternatively equations ( 2.12| ) and (2.f3) imply the existence of a complex-valued function 8 and 
coordinate system (w, z, x, y) —: {w^\xa), as above, such that 

f d d 



dAA' 



^xy^x 



^xx^y 



d 



\dx^ dw^ dx^dx^ dxi 



As a consequence of ( p.f4 ) Q satisfies second heavenly equation 

929 1 a^e 



e,,. 



9^9 



^xx^yy ©xy Or Q^AQj.^ 2 dx^ dx^ DxbOxa 



= 0. 



The dual frame is given by 



dx' 



929 



dx^dxA 



-dw^ 



dw' 



with 9 = defining the flat metric. The Lax pair corresponding to ( ^.22 ) is 

Lq = dy — A(9u) — Qxydy + Qyydx), 

Li = dx + \{dz + Qxxdy - Qxydx)- 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



Both heavenly equations were originally derived by Plebanski p5| from the formulation (2.16). The 
closure condition is used, via Darboux's theorem, to introduce ZJ^, canonical coordinates on the spin 
bundle, holomorphic around A = such that the two-form (2.15) is S(A) = dhUJ^ A dhUOA- The various 
forms of the heavenly equations can be obtained by adapting different coordinates and gauges to these 
forms. 



3 The recursion operator 



I n §§|3.1| the recursion operator R for the anti-self-dual Einstein vacuum equations is constructed. In 
§§3.2 then show that the generating function for Wcj) is automatically a twistor function, and is in fact 
a Cech representative for (p. It is shown that R acts on such a twistor function by multiplication. A 
similar application to the coord inat es used in the heavenly equations yields the coordinate description of 
the twistor space starting. In §§3_^ we show how that the action of the recursio n o perator on space-time 
corresponds to multiplication of the corresponding twistor functions by A. In §§3^ the algebra of hidden 
symmetries of the second he avenly equation is constructed by applying the recursion operator to the 
explicit symmetries. In §§3_^, R is used to build a higher valence Killing spinors corresponding to hidden 
symmetries. In the last subsections examples of the use of the recursion operator are given. 



3.1 The recursion relations 

The recursion operator i? is a raap from the space of hnearised perturbations of the ASDVE equations to 
itself. This can be used to construct the ASDVE hierarchy whose higher flows are generated by acting on 
one of the coordinate flows with the recursion operator R. 

We will identify the space of linearised perturbations to the ASDVE equations with solutions to the 
background coupled wave equations in two ways as follows. 

Lemma 3.1 Let Dn and Dq denote wave operators on the ASD background determined by and & 



respectively. Linearised solutions to ( 2.18| ) and ( 2.22 ) satisfy 



DnSn = 0, De^e = 0. (3.25) 
Proof. In both cases Dg — Wav'^^q' since 

but dag"^'' = for both heavenly coordinate systems. For the flrst equation {dd{fl + (50))^ = v implies 

= {ddn A dd)sn = d{ddn Aid- d)sn) = d * dsn. 

Here * is the Hodge star operator corresponding to g. For the second equation we make use of the tetrad 



(|2.21| ) and perform coordinate calculations. 

□ 



From now on we identify tangent spaces to the spaces of solutions to (2.18) and (2.22) with the space 
of solutions to the curved background wave equation, Wg. We will define the recursion operator on the 
space Wg. 

The above lemma shows that we can consider a linearised perturbation as an element of Wg in two 
ways. These two will be related by the square of the recursion operator. The linearised vacuum metrics 
corresponding to dfl and SO are 

h^AA'BB' — I'iA'OB')'^ (Al'"^ B)0'^^^ h^^AA'BB' — OA'Ob'^ AO'^ BO'SO- 

where o"^' = (1,0) and t-^' = (0, 1) are the constant spin frame associated to the null tetrads given above. 
Given (j) £ Wg we use the first of these equations to find . If we put the perturbation obtained in 
this way on the LHS of the second equation and add an appropriate gauge term we obtain (j)' - the new 
element of Wg that provides the 6Q which gives rise to 

Kl = hi, + W^aVh). (3.26) 

To extract the recursion relations we must find V such that h/ aa'bb' — ^ (AA'Vbb') = oa'Ob'Xab- Take 
Vbb' = OB'V_Bi'(5ri, which gives 

^{AA'VbB') = -'-(A'Os')V(A0'VB)i/(5fi + OA'Ob'V AV^ Bl'Sfl. 



This reduces (B.26) to 

V^i'Viji'-Z-^ Vao'Vso'0'- (3.27) 
Definition 3.2 Define the recursion operator R : Wg — > Wg by 

i'^'VAA'^^0^'VAA'R4>, (3.28) 

so formally R — (V^o')^^ ° V^ii' [no summation over the index A). 
Remarks: 



From ( 3.28 ) and from ( 2.11 ) it follows that if 4> belongs to Wg then so does R4>. 



• If R^Sn — 6Q then Sft and 6Q correspond to the same variation in the metric up to gauge. 

• The operator (j) V ao' is over-determined, and its consistency follows from the wave equation on 

• This definition is formal in that in order to invert the operator cf) i—f \/ao'4' we need to specify 
boundary conditions. 



To summarize: 



Proposition 3.3 Let he the space of solutions of the wave equation on the curved ASD background 
given by g. 

(i) Elements ofWg can be identified with linearised perturbations of the heavenly equations. 



(ii) There exists a (formal) map R : Wg — > Wg given by (3.2 



The recursion operator can be generalised to act on solutions to the higher helicity Zero Rest-Mass 
equations on the ASD vacuum backgrounds [|lO) by using Herz potentials. We restrict ourselves to the 
gauge invariant case of left-handed neutrino field ipA on a heavenly background. First note that any 
solution of 

V-^^Va = 

must be of the form Vao"/* where (j) £ Wg. Define the recursion relations 

TZ2Pa-=^ao'R^- (3.29) 

It is easy to see that TZ maps solutions into solutions, although again the definition is formal in that 
boundary conditions are required to eliminate the ambiguities. A conjugate recursion operator TZ will 
play a role in the Hamiltonian formulation in Section ^. 



3.2 The recursion operator and twistor functions 

A twistor function / can be pulled back to the correspondence space F. A function f on F descends to 
twistor space iff L^/ = 0. 

Given cj) G Wg, define, for i e Z, a hierarchy of linear fields, = -RVo- Put '5 J2°^oo't'i^^ ^^'^ 
observe that the recursion equations are equivalent to La'^ — 0. Thus ^E* is a function on the twistor 
space VT. Conversely every solution of La'^' — defined on a neighbourhood of |A| = 1 can be expanded 
in a Laurent series in A with the coefficients forming a series of elements of Wg related by the recursion 
operator. The function 5*, when multiplied by l/(7ro'7ri'), is a Cech representative of the element of 
H^{VT,0{—2)) that corresponds to the solution of the wave equation <j) under the Penrose transform 
(i.e. by integration around |A| — 1). The ambiguity in the inversion of V^o' means that there are many 
such functions 5* that can be obtained from a given (f). However, they are all equivalent as cohomology 
classes. 

It is clear that a series corresponding to i?0 is the function A^^^*. As noted before, R is not completely 
well defined when acting on Wg because of the ambiguity in the inversion of V ao' ■ However, the definition 
R^ = is well defined as a twistor function on VT, but the problem resurfaces when one attempts to 
treat ^'(A) as a representative of a cohomology class since pure gauge elements of the first sheaf cohomology 
group H^{'PT,0{—2)) are mapped to functions defining a non-trivial element of the cohomology. Note, 
however, that with the definition R^i = ^'/A, the action of R is well defined on twistor functions and can 
be iterated without ambiguity. 

We can in this way build coordinate charts on twistor space from those on space-time arising from 
the choices in the Plebanski reductions. Put luq = = {w,z); the surfaces of constant cu^ are twistor 
surfaces. We have that \/^q'UJq = so that in particular V ai'^^o'^o = and if we define Lof = R^uJq 
then we can choose Luf- — for negative i. We define 

oo 

u:^ = Y,u:t\\ (3.30) 

1=0 

We can similarly define Cj^ by Cjq — and choose uof — for i > 0. Note that w'^ and uj^ are solutions 
of La holomorphic around A = and A oo respectively and they can be chosen so that they extend to a 
neighbourhood of the unit disc and a neighbourhood of the complement of the unit disc and can therefore 
be used to provide a patching description of the twistor space. 



3.3 The Penrose transform of hnearised deformations and the recursion op- 
erator 

The recursion operator acts on linearised perturbations of the ASDVE equations. Under the twistor 
correspondence, these correspond to linearised holomorphic deformations of (part of) VT . 

Cover VT by two sets, U and U with |A| < 1 + e on ?7 and |A| > 1 — e on [/ with {uj^ , A) coordinates 
on U and {lu^,\~^) on U. The twistor space VT is then determined by the transition function lu^ = 
ui^iuj^ , TTA') on [/ n [/ which preserves the fibrewise 2-form, Auj^ A da;A|A=const. = do)"^ A d(I'A|A=const.- 



Infinitesimal deformations are given by elements of H^{VT, &), where denotes a sheaf of germs of 
holomorphic vector fields. Let 



defined on the overlap U HU and define a class in H^{VT, 0) that preserves the fibration VT t-^ CP^. 
The corresponding infinitesimal deformation is given by 

w^(t^^,7r^,,t) = (l+ty)(tD^) + 0(i2). (3.31) 

From the globality of S(A) = dcj^ AduJA it follows that y is a Hamiltonian vector field with a Hamiltonian 
/ e H^{'PT, 0{2)) with respect to the symplectic structure E. A finite deformation is given by integrating 

d^^ ,BA df 



= e 



dt duj^ ■ 

from t = to 1. Infinitesimally we can put 

OLOA 

If the ASD metric is determined by Q and then e^^dSf /duj^ , (or more simply 6f) is a linearised defor- 
mation corresponding to SQ £ Wg . 

The recursion operator acts on linearised deformations as follows 



Proposition 3.4 Let R be the recursion operator defined by (3.28). Its twistor counterpart is the multi- 
plication operator 

R5f^^Sf = X-'Sf. (3.33) 

TTO' 

[Note that R acts on 5f without ambiguity; the ambiguity in boundary condition for the definition of R 
on space-time is absorbed into the choice of explicit representative for the cohomology class determined 
by Sf.] 

Proof. Pull back 6f to the primed spin bundle on which it is a coboundary so that 

6fi7:A',x'')^hinA',x'')-h{nA',x'') (3.34) 

where h and h are holomorphic on U and U respectively (here we abuse notation and denote by U and 
U the open sets on the spin bundle that are the preimage of U and U on twistor space) . A choice for the 
splitting ( 3.34 ) is given by 

^ ^ :5f{PE')PD'dp''\ (3.35) 



2ni Jy {p^' 'Kc'){p^' ob'Y 

1 / (tT^'oaO^ .r, N . D' 



27ri Jf {p^' 'nc'){p^' ob'Y 

Here pA' are homogeneous coordinates of CP^ pulled back to the spin bundle. The contours F and F are 
homologous to the equator of CP^ in C7 n C/ and are such that F — F surrounds the point pA' — t^A' ■ 

The functions h and h are homogeneous of degree 1 in -ka' and do not descend to "PT, whereas their 
difference does so that 

TT^'VAA'h = TT^'^AA'h = Tr^'7T'''7T^'^AA'B'C' (3.36) 

where the first equality shows that the LHS is global with homogeneity degree 2 and implies the second 
equality for some "^aa'D'c which will be the third potential for a linearised ASD Weyl spinor. Tiaa'B'C 
is in general defined modulo terms of the form V a{A'1b'C') but this gauge freedom is partially fixed by 
choosing the integral representation above; h vanishes to third order at tta' = oa' and direct differentiation, 
using Vaa'S/ = Pa'S/a for some S/a, gives T,aa'B'C' = oa'Ob'Oc'"^ aq'SQ where 



S^^^f , .. Po'dp"" . (3.37) 

27r« Jr {p^ OB') 

This is consistent with the Plebanski gauge choices (there is also a gauge freedom in 5Q arising from coho- 
mology freedom in 5f which we shall describe in the next subsection.) The condition \7 AiD'"^^ A' b'C) — 
follows from equation ( 3.36 ) which, with the Plebanski gauge choice, implies SO £ Wg. Thus we obtain 
a twistor integral formula for the linearisation of the second heavenly equation. 



Now recall formula ( 3.28 ) defining R. Let RSf be the twistor function corresponding to R6Q by ( 3.37 ). 
The recursion relations yield 



RSf A , D' I SfA , D' 



so RSf = \-^5f. 

□ 

Let 5VL be the linearisation of the first heavenly potential. From R^Sft — 6Q it follows that 

1 /" S f 

'''-^^iiP.o-'npB'.-y''"'p"'■ 

3.4 Hidden symmetry algebra 

The ASDVE equations in the Plebahski forms have a residual coordinate symmetry. This consists of 
area preserving diffeomorphisms in the coordinates together with some extra transformations that 
depend on whether one is reducing to the first or second form. By regarding the infinitesimal forms of 
these transformations as linearised perturbations and acting on them using the recursion operator, the 
coordinate (passive) symmetries can be extended to give 'hidden' (active) symmetries of the heavenly 
equations. Formulae ( |3.37 ) and ( [3.33| ) can be used to recover the known relations (see for example [ p8| ) 



of the hidden symmetry algebra of the heavenly equations. We deal with the second equation as the case 
of the first equation was investigated by other methods ||2^ . 

Let M be a volume preserving vector field on A4. Define (5"^Vaa' ■— [M,W aa']- This is a pure 
gauge transformation corresponding to addition of C^g to the space-time metric and preserves the field 
equations. Note that 

['>Ii^S%]Vaa' := ^"m,jv]'^aa'- 

Once a Plebahski coordinate system and reduced equations have been obtained, the reduced equation will 
not be invariant under all the SDiff(Al) transformations. The second form will be preserved if we restrict 
ourselves to transformations which preserve the SD two- forms ^ = dwA Adw"^ and T,^ ^ — dxA /\dw^ . 
The conditions Cm^°'°' = Cm'S°'^' = imply that M is given by 

dh d f do R d'^h \ d 
M=- -— r + Ur^-a;^ 



dwAdw^ \dwA OwaOw^ J dx^ 

where h — h{w^) and g — g{'w'^). The space-time is now viewed as a cotangent bundle A4 = T*N'^ with 
being coordinates on a two-dimensional complex manifold A/"^. The full SDiff(A^) symmetry breaks 
down to the semi-direct product of SDiff(A/'^), which acts on by a Lie lift, with r(A/'^, O) which acts on 
M. by translations of the zero section by the exterior derivatives of functions on N'^ . Let i5m6 correspond 
to 6%j\Iaa' by 

Om^ Al 

The 'pure gauge' elements are 



dx^dx^ dxB 



5\j<d = F + xaG^ + xaXb^ ^ VxaxbXc 



dwAdwB dwAdwBdwc 
dg 99 ^ dh de g d^h 99 



^"ttT^TJZT (3-38) 



dwA dx^ dwA dw^ dwAdw^ dx^ 

where F, G^, h are functions of only. 

The above symmetries can be seen to arise from symmetries on twistor space as follows. Since we 
have the symplectic form S — dw^ A dtOA on the fibres of p : VT — > CP^ , a symmetry is a holomorphic 
diffeomorphism of the set U that restricts to a canonical transformation on each fibre. Let H — H{x°', A) = 
Si^o Hamiltonian for an infinitesimal such transformation pulled back to the projective 

spin bundle. The functions hi depend o n spa ce time coordinates only. In particular ho and hi give 



h and g from the previous construction ( 3.38 ). This can be seen by calculating how 9 transforms if 
u!^ = w"^ + Xx^ + X^dQ/dxA + ... — > <^^- Now 9 is treated as an object on the first jet bundle of a 
fixed fibre of VT and it determines the structure of the second jet. 

These symmetries take a solution to an equivalent solution. The recursion operator can be used to 
define an algebra of 'hidden symmetries' that take one solution to a different one as follows. 



Let SifQ be an expression of the form ( 3.38 ) which also satisfies DgJ^jG = 0. We set 



Proposition 3.5 Generators of the hidden symmetry algebra of the second heavenly equation satisfy the 
relation 

[6m\Sn']^Simm]'^'- (3.39) 

Proof. This can be proved directly by showing that the ambiguities in R can be chosen so that Ro5m — 
Sm o R- It is perhaps more informative to prove it by its action on twistor functions. 

Let S\.jf be the twistor function corresponding to S\jQ (by ( |3.37| )) treated as an element of T{U f] 
U,0{2)) rather than H^(VT,0{2)). Define [5\i,5%] by 



where the Poisson bracket is calculated with respect to a canonical Poisson structure on VT. From 



Proposition (|3.33|) it follows that 



2m J (vro'j^ 

as required. 

□ 

3.5 Recursion procedure for Killing spinors 

Let {Ai,g) be an ASD vacuum space. We say that is a Killing spinor of type (0,n) if 

V^iA'LB',...B'j - 0. (3.40) 
Killing spinors of type (0, n) give rise to Killing spinors of type (1, ti — 1) by 

'^^A'Lb[...B!,^ = £a'(b[K^ B'^...B'J- 

In an ASD vacuum, K^^'^- -^" is also a Killing spinor 
Put (for i — 0, n) 

L, := t^i...i^-o^'+i...o^"Ls;...S' , 



and contract ( |3.40[ ) with l^k..l^''o^'^+^ ...o^'"-+^ to obtain 

iVAvLi-i = -(n - z + l)VAQ'Li, z = 0, n - 1. 



We make use of the recursion relations ( 3.28| ): 

— i 

— — :R[Li^i) = Li. 

n + 1 — z 

This leads to a general formula for Killing spinors (with V^o'^o = 0) 



L^^{-lr\.] R\Lo), LB'^B^...B'^^J2''iB^■■■0B'^B'^,■■■^B')L^ (3.41) 



and equation ( 3.4C ) is then satisfied iff i? ^Lq = = 0. 



3.6 Example 1 



Let us demonstrate how to use the recursion procedure to find metrics with hidden symmetries. Let 
dt„^ := (j)n be a Unearisation of the first heavenly equation. We have R : z — > ft^ = dt^i^. Look 
for solutions to ( 2.1^ ) with an additional constraint St^fi = 0. The recursion relations (3.28) imply 
rtwz — ^ww = 0, therefore 

il{w, z, w, z) ~ wq{w, z) + P{z, w, z). 
The heavenly equation yields dq A dP A dz = dz A dw A dz. With the definition dzP = p the metric is 

ds^ = 2dwdq + 2dzdp + /dz^ 

where / = —2Pzz- We adopt {w,z,q,p) as a new coordinate system. Heavenly equations imply that 
/ ~ f{q, z) is an arbitrary function of two variables. These are the null ASD plane wave solutions. 



3.7 Example 2 



Now we shall illustrate the Propositions p.3| and 3^ with the example of the Sparling-Tod solution [£7| 
The coordinate formulae for the pull back of twistor functions are: 



Consider 



- Xy- X^e^ + X^Qz 

Xx - A^e,, - A^e,,, 



e 



wx + zy 



(3.42) 



(3.43) 



where cr = const. It satisfies both the line ar and the nonlinear part of ( 2.22 ). 

The fiat case: First we shall treat ( 3.43 ), with fi = 1, as a solution 0o to the wave equation on the flat 
background. The recursion relations are 



y 



{wx + zyY ' 



— X 



(wx + zyY 



They have a solution (pi :— i?0o — {~'y/w)(j>o. More generally we find that 



wx 



zy 



(3.44) 



The last formula can be also found using twistor methods. The twis tor function corresponding to (f>o is 

the twistor function corresponding 



Ay and fii = z — Xx. By Proposition 3.33 



l/(/i°/i"'^), where fig 

to (fin is A~"/(/i°/x^). This can be seen by applying the formula (3.37) and computing the residue at the 
pole A = —w/y. It is interesting to ask whether any 0„ (apart from (/jq) is a solution to the heavenly 
equation. Inserting 8 = 0„ to (2.22) yields n = or n = 2. We parenthetically mention that 02 yields 



(by formula ( ^.23| )) a metric of type D which is conformal to the Eguchi-Hanson solution. 
The curved case. Now let Q given by (3.43) determine the curved metric 



ds^ — 2du'da; + 2dzdy + 4:(j{wx + zy) '^[wdz — zdwY 



(3.45) 



The recursion relations 



dy{R(t)) = [dw - Qxydy 



Qyydx) 



'dx{R(j)) = [dz + Qxxdy - <dxydx)<l> 



are 



'dx{Rip) = {dz + 2aw{wx + zy) {wdx - zdy))^) , 



dy{RiP) 



{dw + 2az{wx + zy) ^ {wdx — zdy))ip , 



where ip satisfies 



□eV" — "^{dxdw + dydz + 2a{wx + zy) ^{z^dx^ + w'^dy^ — 2wzdxdy))ip — 



..2a 2 



(3.46) 



One solution to the last equation is tpi — {wx + zy) ^. We apply the recursion relations to find the 
sequence of linearised solutions 



4>2 



V w/ wx 4 " " 

n 



a 



wlwx + zy' 3 {wx + zy)^ 

k 

k — n 



1 



y 



w/ wx + zy 



k=0 



To find A^s note that the recursion relations imply 



R 



— — ^ (wx + zyY 



1 - 



-o- {wx + zy) ——7 (wx + zy) 

w/ \ wJ J + 2/ \ w 

This yields a recursive formula 

which determines the algebraic (as opposed to the differential) recursion relations between V'n and 4'n+i- 



It can be checked that functions ipn indeed satisfy (3.46). Notice that if ct = (flat background) then we 



recover (3.44). We can also find the inhomogeneous twistor coordinates pulled back to T 

oc n ^ 

OO Tl ^ 



n=0 k=0 



where 



- - 2a— — = 1, i?i = 0, B.-i = 0, = 0...n. 



fc+ 1 

The polynomials solve L^(/i^) = 0, where now 

Lo = — A3tu — 2A<Tz^(w2; + zj/)^'^9a; + (1 + 2AfTwz(wz + z2/)~'^)9y, 
Li = \dz + [1 ~ 2X<Twz{wz + zyy^)dx + 2Xcrw'^{wz + zy)~^)dy. 

4 Hierarchies for the ASD vacuum equations 

The hidden symmetries corresponding to higher flows associated to translations along the coordinate 
vector fields give 'higher flows' of a hierarchy. This yields a hierarchy of flows of the anti-self-dual Einstein 
vacuum equations. We first give this for the equations in their second heavenly form but then give the 
equations in the form of consistency conditions for a Lax system of vector fields generalizing equations 



2.11. The nonlinear graviton construction generalizes to give a construction for the corresponding system 



of equations and is presented in §§4.2. In §§|4.3| the geometric structure of solutions to the truncated 



hierarchy are explored in further detail. Finally in §§4.4 infinitesimal deformations are studied 



4.1 Hierarchies for the heavenly equations 

The generators of higher flows are first obtained by applying powers of the recursion operator to the 
linearised perturbations corresponding to the evolution along coordinate vector fields. This embeds the 
second heavenly equation into an infinite system of over-determined, but consistent, PDEs (which we will 
truncate at some arbitrary but finite level). These equations in turn can be naturally embedded into 
a system of equations that are the consistency conditions for an associated linear system that extends 



( 2.11 ). We shall discuss here the hierarchy for the second Plebahski form; that for the first arises from a 
different coordinate and gauge choice. 

Introduce the coordinates x^"^, where for i = 0, 1, x^^ = x^^ are the original coordinates on A4, and 
for 1 < I < n,x^^ are the parameters for the new flows (with 2n — 2 dimensional parameter space X). 
The propagation of O along these parameters is determined by the recursion relations 

dy{dBi+iO) = [dio - Oxydy + Qyydx)dBiQ , 
-dx{dBi+i&) = [dz + Qxxdy ~ Sxydx)dBiQ , 
or ^Ao{^B^+le) = {^Al+^co^AoQ^^o)^B^e■ (4.48) 

However, we will take the hierarchy to be the system (containing the above when j = 1) 

dAidBj-iQ - dB^dAi^iQ + {dAi-iQ,dBj-iQ}yx ^ 0, i,j^l...n. (4.49) 
Here {..., ■■■}yx is the Poisson bracket with respect to the Poisson structure d/dx^ A d/dxA — 2dx A dy. 



Lemma 4.1 The linear system for equations (4.49) is 

LaiS = {-XDAi+1 + SAi)s = 0, z = 0, n - 1, 

where 

1. s := s(x^% A) is a function on a spin bundle (a -bundle) over = A1 x X, 

2. DAi+i 9^1+1 + [dAi, V], {V — e^^ dAoQdBo) o.^d Sai '■— Oai are 4n vector fields on M . 



(4.50) 



Proof. This follows by direct calculation. The compatibility conditions for (4.50) are: 

[Dm+uDb-j+i] = Q, (4.51) 

[<5A»,fej]=0, (4.52) 

[Dm+u 5Bj\ - [Dbj+u 5m] = 0. (4.53) 

It is straightforward to see that equations ( 4.52| ) and (4.53) hold identically with the above definitions 
and (4.51) is equivalent to (4.49). 

□ 



As a converse to this lemma, we will see in §§T2 using the twistor correspondence, that given the Lax 
system above, in which the vector fields Dai and Saj are volume preserving vector fields, then coordinate 
and gauge choices can be made so that the Lax system takes on the above form. 

4.1.1 Spinor notation 

The above can also be represented in a spinorial formulation that will be useful later. We introduce the 
spinor indexed coordinates x^^'^^^^^" = x'^^'^i - '^"' on M which correspond to the x^^ by 

\™AAi A' ...A' 



'OAi---OA'/A;^i---iA^(-l)" 

The vector fields DAi+i and Sai are then represented by the An vector fields on N, DAA[{Ai,...A'^) where 

DAA[i = ''^'^■■■I^^^O^'-+^...0^'^Daai^(ai^ DaiH = DAi+l, Daq'i = SAi 

and La(a'^...a'J = t^"^'^ Daa[(A'^...A'J, Lai = tt'^'^Daa'^i- In the adopted gauge 

Dao'A'^...A'^ — dAO' A'^...A'^i Dai'A'^...A'^ — dAl' A'^...A'^^ + [Qao' A'^...A'^iy]- 

In what follows we will often be interested in a(A{A'^...A' )j the symmetric part of Daa{A'^...A'^- 



■At 



D 



A{A{A'^...A'Jl- 



(4.54) 
(4.55) 



1 i 
= -{iDAi'i-i + [n ~ i)DAa'i) = 9^,; + -[dAi-i,V]. 
n n 

Put Dao'...o' = Qa- The 2n + 2 vector fields 

^ AA{...A'^ = {Qa^^ AO'V A'^...A'^_^,DAn} 

span T*N. 

4.2 The twistor space for the hierarchy 

The twistor space VT for a solution to the hierarchy associated to the Lax system on M as above is 
obtained by factoring the spin bundle A/'x CP"'^ by the twistor distribution (Lax system) Lai- This clearly 
has a projection q : Af x CP^ i-^ VT and we have a double fibration 



TV X CP^ 



v/ 



VT 



Since the twistor distribution is tangent to the fibres of A/" x CP^ ^ CP^ , twistor space inherits the 
projection /x : VT ^ CP"'^. The twistor space for the hierarchy is three-dimensional as for the ordinary 
hyper-Kahler equations, but has a different topology. We have 



Lemma 4.2 The holomorphic curves g(CP^) where CP^ = p ^x, x G Af, have normal bundle N = 
0{n)®0{n). 

Proof. To see this, note that N can be identified with the quotient p* {T^M) / {spanLAi} , i = 1, . . . ,n. 
In their homogeneous form the operators Lai have weight 1, so the distribution spanned by them is 
isomorphic to the bundle C^" ® 0{—l). The definition of the normal bundle as a quotient gives 

^ C^" ® 0{-l) ^ C2"+2 ^ N 

and we see, by taking determinants that the image is 0{n + a) 0{n — a) for some a. We see that a — Q 
as the last map, in the spinor notation introduced at the end of the last section, is given explicitly by 
yAA\...A'^ l/^^'i---^"7r^/^ . . . -KA'^ clearly projecting onto 0{n) ® 0{n). 

□ 

A final structure that VT possesses is a skew form S taking values in 0{2n) on the fibres of the 
projection /i. This arises from the fact that the vector fields of the distribution preserve the coordinate 
volume form v on M \n the given coordinates system. Furthermore, the Lax system commutes exactly 
[La/, Lsj] = so that 

E — v{-, •, LqI, • ■ ■ ■ -^On, ill, • • ■ , Lin) 

descends to the fibres of VT i~> CP^ and clearly has weight 2n as each of the Lai has weight one. 

Thus we see that, given a solution to the hyperkahler hierarchy in the form of a commuting Lax 
system, we can produce a twistor space with the above structures. Now we shall prove the main result 
of this section and demonstrate that, given VT , with the above structures, we can construct M ( as the 



moduli space of rational curve s in V T) which is naturally equipped with a function O satisfying ( 4.4£ ) 
and with the Lax distribution ( 4.5C| ). 



Proposition 4.3 Let VT be a 3 dimensional complex manifold with the following structures 

1) a projection p : VT — > CF^ , 

2) a section s : CP"' ^ VT of p with normal bundle 0{n) © 0{n), 

3) a non-degenerate 2-form S on the fibres of pL, with values in the pullback from CP"'^ of 0{2n). 

Let Af be the moduli space of sections that are deformations of the section s given in (2). Then M is 
2n + 2 dimensional and 

a) There exists c oordi nates, x^^ , A — 0,1, and i — 0, . . . ,n and a function Q : JV — > C on Af such 
that equation ( 4.49| ) is satisfied. 



b ) The moduli space Af of sections is equipped with 

— a factorisation of the tangent bundle TAf — (E) Q"S^ , 

— a 2n- dimensional distribution on the 'spin bundle' D C T{Af x CP^) that is tangent to the fibres 
of r over and, as a bundle on Af x CP^ has an identification with 0{—l) (8) Saa'...a' so 



that the linear system can be written as in equation ( 4.50 ). 

This correspondence is stable under small perturbations of the complex structure on VT preserving (1) 
and (3). 

Proof: The first claim, that Af has dimension 2n + 2 follows from Kodaira theory as dim iJ" (CP^ , N) = 
2n + 2 and dimiJi(CP\ iV) ^ dimi7i(CP\ EndTV) = 0. 

Proof of (a): wc first start by defining homogeneous coordinates on VT. These are coordinates on T, 
the total space of the pullback from CP^ of the tautological line bundle 0{—l). Let -ka' be homogeneous 
coordinates on CP^ pulled back to T and let oj^ be local coordinates on T chosen on a neighbourhood of 
M^^{^o' = 0} that are homogeneous of degree n and canonical so that E = e^sdu;'^ A dw^. We also use 
A ~ ttq' /ttv as an affine coordinate on CP^. Let Lp be the line in VT that corresponds to p € Af and let 
Z e VT lie on Lp. We denote by T the correspondence space VT x Af\zeLp = Af x CP^. (See figure 1 
for the double fibration picture.) 

Pull back the twistor coordinates to and define 2{n + 1) coordinates on Af by 



Figure 1: Double fibration. 

dimF=2n+3 



dimPr=3 




dimAr=2n+2 



where the derivative is along the fibres of !F over J\f. This can alternatively be expressed in afRne 
coordinates on CP^ by expanding the coordinates lu^ pulled back to J-' in powers of A = ttq' / ny : 



(4.56) 



i=0 



where the sf are functions of x"^"*! - "*" and will be useful later. 

The symplectic 2-form E on the fibres of ^, when pulled back to the spin bundle, has expansion in 
powers of A that truncates at order 2n + 1 by globality and homogeneity, so that 

for some symmetric spinor indexed 2-form E^i -'^n^i -^ii . We have 

S(A) A E(A) = 0, dftS(A) = 0. (4.57) 

where in the exterior derivative d^, A is understood to be held constant. 

If we express the forms in terms of the x^"^ and the sf, the closure condition is satisfied identically, 
whereas the truncation condition will give rise to equations on the sf allowing one to express them in 
terms of a function 8(0;"*"* ■■■^n) and to field equations on Q as follows. 

To deduce the existence of <d{x^^i - "^") observe that the vanishing of the coefficient of A^"+^ in 
duj"^ A duJA gives 

^ dsAi A da;^* = d ^ SAida;^* = =^ SAi 

i=0 i=0 



dx 



Ai ■ 



The equations of the hierarchy arise from the vanishing of the coefficient of A^"+^ 

n 

da;-^' A ds'+^ + ds^° A d4 = . 



This leads to the equations (4.49) on O for i, j < n ~ \ 



929 



929 



-CD 



dx^'+WxBj dx^^dxBi+'L 
and further equations that determine s"*""*"^. 



dxCOdx^' dxDOdxBi 



= 



Proof of b). The isomorphism Tf\f ~ ® Q^S^ follows simply from the structure of the normal 
bundle. From Kodaira theory, since the appropriate obstruction groups vanish, we have 

T^AA = r(CP;^,iV^) = S'^®0"S"^' (4.58) 

where iV^; is the normal bundle to the rational curve CP;^ in VT corresponding to the point x G M. The 
bundle on space-time is the Ward transform of 0{—n) ® TyVT where the subscript V denotes the 
sub-bundle of the tangent bundle consisting of vectors up the fibres of /i, the projection to CP^, so that 
= r(CP^, 0(-n) ® TvVT). The bundle S"^' = r(CP\ 0(1)) is canonically trivial. 



Let AA'^---A' — A{A'^---A' ) be the indexed vector field that establishes the isomorphism (4.58) and 
let e-^^'i-^" = e-^^-^'i-O &VL^ ^S"^® ©"S"^' be the dual (inverse) map. 

We now wish to derive the form of the linear system, equations ( 4.50 ). For each fixed tta' = (A, 1) G 
CP"'^ we have a copy of a space-time J\f\. The horizontal (i.e. holding A constant) subspace of T(^x,\){-^ x 
CP^) is spanned by yA{A'...A' )• An element of the normal bundle to the corresponding line CPiJ, consists 
of a a horizontal tangent vector at [x, A) modulo the twistor distribution. Therefore we have the sequence 
of sheaves over CP^ 

— > Dx — > T^M ^ S^® ©(n) — > , 

where is the twistor distribution at x and the map T^J^ — > © 0{n) is given by the contraction 
of elements of T^Af with := e^^'-'"^"TrA[---TTA'^ since annihilates all LsiS in D. Consider the dual 
sequence tensored with 0{—l) to obtain 

OAi^n - 1) ^ T;AA(-1) ^ D:{-1) 0. (4.59) 

From here we would like to extract the Lax distribution 

Laa',...A'^ = t^'^''Daa[A',...A'„ e Saa',...A'^ <^ 0{1) (g> D. 



This can be achieved by globalising (4.5£) in tt'^ . The corresponding long exact sequence of cohomology 
groups yields 

nOAi-n - 1)) r(rW(-l)) riD*i-l)) H\OA{-n - 1)) 

H\T*M{-1)) ... 

which (because T*J\f is a trivial bundle so that 0{—l) T*J\f has no sections or cohomology) reduces to 

— > r{D*{-l)) H\OA{-n - 1)) — > 0. 

From Serre duality we conclude, since D has rank 2n, that the connecting map S is an isomorphism 
6 : T{D*{-1)) Saa',...A'„- Therefore 

SeTiD(g,Oil)(g>SAA'..A') (4.60) 



is a canonically defined object annihilating lu^ given by (4. 56). 
In index notation we can put 

S = Laa'^.-.a'^ = t^^^Daa[A'2...A'^, 

where Laa'2...A'^ = La(a!,...A' )i the second identity follows from the globality of Laa'2...A'^ and the 
Daa[A'2...A'^ are vector fields on M lifted to A/" x CP^ using the product structure. 
It follows from Laa'2...A' = that if tt^ — then -Daom^.../!' a;^" = so 

,bbL..b' d 



■Dao'A^...A;. - '^A0'~A'^..:A'^g^B0'B'2...B'^ ' 

for some matrix -^Ao'li'^'A' ■ This matrix must be invertible by dimension counting. By multiplying 
Laa'.^...A'^ by the inverse of this matrix, we find we can put 

.bb^...b; _ b Bi^ b; 

^AO'A'^.-.A'^ ^ ^A^A'^-'-^A'^- 

Therefore we can take Laa'^.-.A'^ = dAo'A'^.-.A'^ ~ ^Dava'^.-.A'^- Equating the (n — i + l)th and {n + l)th 
powers of A in LaiOJ^ = to zero yields 

DaVA'.^...A'^ — + [dAO' A'.^...A'^iy] 

where V = e abOQ / dx aqO / Oxeq. So finally Laa'.^...A'^ is of the form LAi — dAi — ^{dAi+i + [dAi, V]). 

□ 



4.3 Geometric structures 

If one considers Af = Ai x K as being foliated by four dimensional slices t"^^ = const then structures 
(l)-(3) on PT can be used to define anti-self-dual vacuum metrics on the leaves of the foliation. Consider 
Q{x^^ ,t) where t = {t^'^,i — 2...n}. For each fixed t the function Q satisfies the second heavenly 
equation. The ASD metric on a corresponding four-dimcnsional slice A/t=to is given by 

d.^ ^ 2e^sd.^U.-^' + 2^^J^d.^^'d.-'. 

This metric can be determined from the structure of the 0{n) 0{n) twistor space as follows. 

Fix the first 2n — 2 parameters in the expansion (4.56) so the normal vector W = W^d / duj^ is given 

by 



where SQ = W"^"^' dO/dx^"^' . The metric is 



giU, W) ^ SABSA'B'U^^'W^^' (4.61) 

where ea'B' is a fixed element of A^S"^ and eab & A^S^ is determined by S; recall that = 
r{Lx,0{-\) (g) TvVT). Thus if u'^,v'^ e 5^, then define eabu^v^ = E(u,w) where u,v are the corre- 
sponding weighted vertical vector fields on VT . 

For n odd TJ\f is equipped with a metric with holonomy SL{2, C). For n even, TJ\f is endowed with 
a skew form. They are both given by 

G{U,W) = eABeA[B[-eA'„B'^U^^''-<W^'''^-'''^^. (4.62) 

These are special examples of the paraconformal structures considered by Bailey and Eastwood ||^. 

4.4 Holomorphic deformations and 0{2n) twistor functions 



We wish to consider holomorphic deformations of VT that preserve conditions (1 — 3) of Proposition 43 
which will therefore correspond to perturbations of the hierarchy. 

Let = G'^{uj^,TTA',t) be the standard patching relation for VT and let e <S) H^iVT,0{n)) 
give the infinitesimal deformation 

The globality of the symplectic structure dcjA A dui^ — dcjA A doj^ implies /'^ = e^^df /duj^ where 
f eH\VT,0{2n)). 

Example: if we deform from the flat model using / = (ttq')*"/^"'^^' then the deformation equations 

imply that Q — lo'^uj^ = cj^cj^ is a global twistor function (up to 0{t^)) which persists to all orders 
as e^^dQ/duj^^df /dui^ — 0. The corresponding deformed paraconformal structure admits a symmetry 
corresponding to the global vector field e^^ dQ / du^d / duu^ on VT . 

To see how such 'Hamiltonians' / correspond to variations in the paraconformal structure (or more 
simply 0), we form an indexed element of H^[VT ,0{—l)), and pull it back to A/" x CP^ where it can be 
split uniquely: 

T^A'^-'-T^A'^ Q^AQ^BQ^C ~ IaBCA'^.-.A'^ = ^ABCA'^.-.A'^ — ^ABCA'^.-.A'^- 

where 

1 / f ABC A'... A 



^ABCA',...A',^ = ^J^ pJa'^ P-^P- 

This gives rise to a global field that is symmetric over its indices: 

CaBCDA'2...A'^D'2...D'^ = LdD'^...D'^^ ABCA'^.-.A'^ 

which is given also directly by the integral 

CaBCDA',...A'„D',...D'„ = i PA',-PA'„PD',...PD'„ Q^AqII Q^C Q^D P ' ^P- 



To see how this corresponds to a variation of 0, we introduce a chain of potentials. Use the non- unique 
sphtting — JF^" — JF^" and define a global object of degree 2n + 1 by 

T -r-2n V ^-B; B'D'C, C' 

i^AA'^...A'^-r ^ ^AA'2...A',^B[...B'„C[...C^D{T^ ^•■•TT " TT ^TT i...7r 

It is easy to see that 

V^^^---^'^'E,AA'2...A'^B[...B'^C[...C"^D[ — 0, 

and '^AA'2...A'^B!^...B'^c'[...C!;^D[ is a potential potentials, related to the field by 

CaBCDA'2...A'^D'2...D'^ — ^ DD'^...D'^^C^ """^B^ '^^AA'2...A'^B{...B'^C[...Ci^D{- 

The chain of potentials is 

SQa[b{...b^c[...C'„d{ = oa'^ob[--.ob'^oc[-.-oc^od[SQ 

^AA'^...A'^B{...B'^C[...C'^D[ = 0b[---0b[^0c[---0C'^0d[V AO' A'^...A'^5Q 

HaBA'^...A'^^B[...B'^^D^, = Ob[---Ob'^Od['^B0'^A0'A'2...A'^SQ 

^ABCA'^...A'^D[ = Vco'Vbo' Vao'A^...A'„^0 

CaBCDA'^...A'^D'^...D'^ — ^ Ca'"^ BO'^ AO' A'^...A'^^ DO' D'^...D'^5Q. 

This can be compared with the corresponding chain for n — 1 . 

5 Hamiltonian and Lagrangian formalisms 

In this Section we shall investigate the Lagrangian and Hamiltonian formulations of the hyper-Kahler 
equations in their 'heavenly' forms. The symplectic form on the space of solutions to heavenly equations 
will be derived, and proven to be compatible with a recursion operator. 

Both the first and second heavenly equations admit Lagrangian formulations, and these can be used 
to derive symplectic structures on the solution spaces, which we denote by S. Here, rather than consider 
the equations as a real system of elliptic or ultra- hyperbolic equations, we complexify and consider the 
equations locally as evolving initial data from a 3-dimensional hyper-surface and it is this space of initial 
data that leads to local solutions on a neighbourhood of such a hyper-surface that is denoted by S and is 
endowed with a (conserved) symplectic form. 

For the first equation we have the Lagrangian density 

Cn = n(v-^{dmf) = (5.63) 



and for the second equation 



ie{e„ Qy}.y-\{Q^Qn,+QyQz)y- (5.64) 



Note that e^^ A e\ can be replaced by da; A dy in the second Lagrangian as it is multiplied by Aw A dz. 

If the field equations are assumed, the variation of these Lagrangians will yield only a boundary term. 
Starting with the first equation, this defines a potential one-form P on the solution space S and hence a 
symplectic structure H — dP on S. Starting with the second we find a symplectic structure with the same 
expression on perturbations 5Q as we had for 5^1. However, since their relation to perturbations of the 
hyper-Kahler structure are different, they define different symplectic structures on S. These are related 
by the recursion operator since we have R^Sft — SO from above. In order to see that these structures yield 
the usual bi-Hamiltonian framework, we will need to show that these symplectic structures are compatible 
with the recursion operator in the sense that S},{R(j), </>') — $1(0, i?(/)'). 

We shall demonstrate this using the first heavenly formulation which is easier as one can use identities 
from Kahler geometry. (The derivation of the symplectic structure from the second Lagrangian will be 
done in coordinates, since the useful relation between the Hodge star and the Kahler structure is missing 
in this case.) 

Proposition 5.1 The symplectic form on the space of solutions S derived from the boundary term in the 
variational principle for the first Lagrangian is 

n{Sin,S2n) = - [ Sin*d{S2n)-S2n*d{5in). (5.65) 

3 JsM 



Proof. Varying ( 5.63|) we obtain 

We use the identities d{d — d) ~ 2dd, uj A Jid = ddil A (9 — 9) = *d and the field equation to obtain 

SL = -\ddn A {snd{d - d)n - nd{d ~ d)6n) 

= \dA{Sn) - \ddn{- * ddn{d - d)6n{d - d)n + *ddn{d - d)n{d - d)sn) 

= \dA{5n) where A{6n) = n*d6n-sn* dn. 

o 

Define the one form on S 

[ A{5Q). 

JsM 

The symplectic structure ft is the (functional) exterior derivative of P 

n{din,S2n) = di{p{S2n)) ^s^iPiSin)) ~ p{[din,S2n]) 

2 



/ 6in*d{d2^) - S2n*d{Sin). □ 

JSM 



Thus O coincides with the symplectic form on the solution space to the wave equation on the ASD 
vacuum background. 

The existence of the recursion operator allows the construction of an infinite sequence of symplectic 
structures. The key property we need is the following 

Proposition 5.2 Let (j), (f)' G Wg and let ft be given by ( 3.65| ). Then 

n{R(j), 4>') = R(f)'). (5.66) 

We first prove a technical lemma: 
Lemma 5.3 The following identities hold 

ui /\d(j) = -a A dR<j), lu A d2<p ^ a A d2R(j), (5.67) 
uj A d2R(f) = -a A dQ(j), uj A dRcj) ^ a A dcj). 

Proof. From the definitions of Y.^ ^ and d^, it follows that 

E^'s' A 9g; = S^'I^' A (5.68) 
(recall that d^, = e^^ ® Qaa' ) which yields 

Lu Ad — a A 82, ujAd = —aAdo, 
uj Ada ~ a Ad, uj A 82 — —a A 82, aAd = aAd — 0. 



Multiplying ( 3.28 ) by combinations of spin co- frame we get an equivalent definition of the recursion 
operator 

d^! <i> = d^! R(j) (5.69) 
which is equivalent to dcj) = d2R4> or = dR(j)- These formulae give the desired result. 



□ 



Proof of Proposition 5.2. The proof uses a (formal) application of Stokes' theorem: 

r2((/>, 0') = ( d(?!)' d0 

JSM 

ijj A {<j)d(j)' - <j)d(j)' - (t)'d(t) + (jj'dct)) = / ujA {(pdcj)' + 4,' dcj)' - 2(l)d(p' - 2(1)' dcj)) 

SM JSM 

-2 [ uuA {<l)d(l)' + (I)' defy) = 2 ( uu A {(fy'dcj) + (l)d<j)'). 

JSM JSM 



From ( 5.7C ) and from ( 5.67 ) we have 



Rcj)') = - [ cuA {(jjdRct)' + Rcji'dct)) 

JSM 

and analogously 



SM 



Aa+ I R(t)'dR4>ha 
'sM 



n{R<j), (j)') = / (f)'d(j) ha- R(t)dR(t)' A a. 

JSM JSM 

Equality ( 5.66| ) is achieved by subtracting the integral of A{4>4>') A a ~ d{R(f>R(j)') A a and applying Stokes' 
theorem. 



□ 



This property guarantees that the bilinear forms 



(5.70) 



are skew. Furthermore they are symplectic and lead to the bi-Hamiltonian formulation. In this context 
formula ( 5.6(]| ) and the closure condition for 17*^ are an algebraic consequence of the fact that R comes 
from two Poisson structures. Using the theory of bi-Hamiltonian systems one can now go on to prove 
that the flows constructed by application of R to some standard flow commute. 

To develop the bi-Hamiltonian theory, we would like to write the heavenly equations in Hamiltonian 
form. However the Legendre transform becomes singular for the coordinate flows associated to the co- 
ordinates we have chosen since they are, at least in the Minkowski space limit, null coordinates. One 
possibility is to develop a Hamiltonian formalism based on such null hyper-surfaces. We shall adopt a 
different approach and reformulate the second heavenly equation as a first order system. 
Define 6 



-Qx and formally rewrite the second heavenly equation (2.22) as 



dn,4> = Tlidy(, 



where TZ ^ {dz + {<j>, ...jyx) ° ^ = Vn' o V 



10' 



(5.71) 



It is therefore a conjugated operator TZ (defined by ( |3.29 )), acting on solutions to the zero-rest-mass 
equations, and plays the role of the recursion operator. Flows of the sub-hierarchy [Li,Loi] = are 



and the Hamiltonian for the first nontrivial fiow is 



Hi 



-dx Ady A dz. 



Higher Hamiltonians iJ„ can in principle be constructed using the operator R. 
developed explicit formulae for these Hn- 



However, we have not 



5.1 A local bi-Hamiltonian form for the hierarchy 

To end this section, we express the equations of the second heavenly hierarchy ( 4.49| ) in a compact 



form, and then write it as a (formal) bi-Hamiltonian system on the spin bundle. This will be a rather 
different framework from that given above in that the Hamiltonian structure will in effect be local to 
the x^^ plane as opposed to a field theoretic formulation — it is the gravitational analogue of that given 
for the Bogomolny equations in except that no symmetries are required here (in effect because ASD 
gravity can be expressed as ASD Yang-Mills with two symmetries but with gauge group the group of area 
preserving diffeomorphisms). This formulation is therefore presented merely as a curiousity. 
Define the jth truncation of ui^ to be 



m—1 

where 9"^* = e^^d/dx^"^. (Note that this is truncated at both ends, although the truncation at the lower 
end and multiplication by a power of A is inessential.) 

Lemma 5.4 The truncated heavenly hierarchy is equivalent to 



(5.72) 



Proof. First observe that one can sum the Lax system to obtain 



= ^^Qaj + {UJA] , ■}yx 



where {/, ■}yx = s'^^dcofdoa- 

Thus, since LAiOJ^ — 0, we have 



which yields the desired answer. 



□ 

For the remainder of this section, we shaU fix the values of the spinor indices to be ^ = and B — 1. Set 

dj :— dij, ^' :— uj'^ , and ^|Jj ujij. 



Equation ( 5.72| ) takes the form 



which we rewrite as 



9,*(A) = {*(A),A-^V,(A)},. 

5,* = 2?^. (5.73) 

Here V :— {^'(A), ■■■}yx — X^i^o ■^"i'^™ A-dependent Poisson structure, Vq = dx and = [9„i_i, V] = 
Dom - dom for m > 0. 
The Hamiltonians are 

6 Outlook - examples with higher symmetries 

This section motivates the study of solutions to heavenly equations which are invariant under some hidden 
symmetries, e.g. along the higher flows. More generally, one can consider solutions to the hyper-Kahler 
equations without symmetries, but whose hierarchies do admit symmetries. 

In a subsequent paper we shall give a general construction of such metrics based on a generalisation of 
||30|] . We consider the case in which the twistor spaces have a globally defined twistor function homogeneous 
of degree n+1. This implies that the metric admits a Killing spinor (some solutions with this property are 
given by §). Global sections Q € H°{C¥\0{n + l)) on non-deformed twistor space fi-.VT — > CP^ will 
be classified and Q-preserving deformations of the complex structure of a neighbourhood of an 0(1)©0(1) 
section of ^ will be studied. The cohomology classes determining the deformation will depend on the fibre 
coordinates of ^ only via Q. The canonical forms of patching functions can be derived to give explicit 
solutions to anti-self-dual ASD vacuum Einstein equation. 

There are also further details of the bi-Hamiltonian structure that could usefully be clarified. 
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